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Abstract 



Linear systems of neutral type are considered using the infinite di- 
mensional approach. The main problems are asymptotic, non-exponential 
. stability, exact controllability and regular asymptotic stabilizability. The 

main tools are the moment problem approach, the Riesz basis of invariant 
subspaces and the Riesz basis of family of exponentials. 

^ '. 1 Introduction 

I Many applied problems from physics, mechanics, biology, and other fields can 

ly^ ■ be described by delay differential equations. A large class of such systems are 

systems of neutral type. In this paper, we consider a general class of neutral 
systems with distributed delays given by the equation 

O; \ -[z{t)-Kzt]^Lzt + Bu{t), t > 0, 

where Zt : [—1,0] C" is the history of z defined by zt{s) = z{t + s). The 



X 

I difference and delay operators K and L, respectively, are defined by 

x/-A_i/(-i) and Lf = l°^A2{e)-^fie)de + J°^A3{e)f{e)de 

for / £ 1, 0], C"), where A_i is a constant nxn matrix, A2, A3 are nxn 

matrices whose elements belong to L2{~1, 0), and i? is a constant nxr matrix. 
A more general case may be when 

Kf= f dfi{e)f{0), /eC([-l,0],C"), 
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where is a matrix- valued function of bounded variation and continuous at zero. 
But we limit ourself to above mentioned case when Kf — 1) because 

that the property of system are mainly characterized by the structure of the 
matrix A^i. For more general forms of the operator K it may expected that 
the situation is analogous, but it is not clear how the properties studied here 
may be connected. Anyway, it is a domain for further investigation. Distributed 
delay may arise in the natural modeling or after some feedback. Our purpose 
is to investigate the problems of asymptotic stability, of stabilizability by linear 
feedback and of exact controllability. 

For that problems, the neutral type systems are less studied that the re- 
tarded systems when K = 0. The difficulties are related to the following par- 
ticular properties of neutral type systems: there may exist an infinite number 
of eigenvalues in the right half plane, in particular near the imaginary axis; the 
choice of the phase-space is crucial, in contrast to the case of retarded func- 
tional differential equations where solutions are more smooth than the initial 
data; some feedback with may change the structure of the system, etc. 



2 The operator model 

In [13] and several other works, the framework is based on the description of neu- 
tral type systems in the space of continuous functions C([— 1, 0]; C"). However, 
for several control problem, the Hilbert space structure is more convenient in 
the study of our class of systems. In Hilbert spaces one can use the fundamental 
tool of Riesz basis (or orthonormal basis modulo a bounded isomorphism). We 
consider the operator model of neutral type systems introduced by Burns and 
al. in product spaces. This approach was also used in [JT] for the construction 
of a spectral model. In 04] the authors consider the particular case of discrete 
delay, which served as a model in [28^ [30] to characterize the stabilizability of a 
class of systems of neutral type. 

The state space is M2(-1,0;C") = C" x i2(-l,0;C"), briefly Ma, and 
permits H]) to be rewritten as 

= Axit) + Bu{t), x{t) = , (2) 

where the operators A and B are deflned by 

A (y(i)\ f f_,A,{9)zt{9)de + f_^A,{9)zt{9)de\ (Bu 

dz,(0)/d^ )' """-(.O 

The domain of A is given by 

V{A) = {{y,z{-)) : z e i^i(-l,0;C"),y = z(0) - A^^-l)} C M2 

and the operator A is the infinitesimal generator of a Co-semigroup e"^*. The 
relation between the solutions of the delay system ([1]) and the system ([2]) is 
zt{9)=z{t + 9). 

In the particular case when A2{9) = A^{9) = 0, we use the notation A for 
A. The properties of A can be expressed mainly in terms of the properties of 
matrix A_i only. Some important properties of A are close to those of A. 
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3 Spectral analysis 



Let us denote by ^i, fii ^ if i ^ j, the eigenvalues of A^i and the 

dimensions of their rootspaces (generahzed eigenspaces) hj pi, ...,p£, J2k=iPk — 
n. Consider the pointsAm"* = In |/i,n| + i(arg/i,„ + 2TTk), m = 1, ..,£;k G Z and 
the circles L^fi} of fixed radius r < ro = imin{|Am'' — '^["'''l, ('ti, fc) ^ (hi)} 
centered at Am''- 

Theorem 3.1. The spectrum of A consists of the eigenvalues only which are 
the roots of the equation dot A(A) = 0, where 

/o M 
e^'A2{s)ds + j e^'A3{s)ds. (4) 

The corresponding eigenvectors of A are = (C — e^'^v4_iC, e'^^C) , with C G 
KerA(A). 

There exists Ni such that for any \k\ > Ni, the total multiplicity of the roots 
of the equation dct A(A) = 0, contained in the circle Lm\ equals Pm- 

The description of the location of the spectrum of we use Rouche theorem. 



3.1 Basis of invariant subspaces 

The most desired situation for concrete systems is to have a Riesz basis formed 
by eigenvectors of A or, at least, by generalized eigenvectors. In more general 
situations, one studies the existence of basises formed by subspaces. We remind 
that a sequence of nonzero subspaces {V^}"" of the space V is called basis (of 
subspaces) of the space V , if any vector x £ V can be uniquely presented as 
X — J^'kLi^k, where Xk G T4, k — 1,2,.. We say that the basis {Vfe}°° is 
orthogonal if Vi is orthogonal to Vj when i ^ j. A basis {Vk} of subspaces is 
called a Riesz basis if there are an orthogonal basis of subspaces {M^fe} and a 
linear bounded invertible operator R, such that RVk = Wk- 

The best "candidates" to form the basis of subspaces are generalized eigen- 
spaces of the generator of a semigroup, but there are simple examples (see 
Example 13.31 below) showing that generalized eigenspaces do not form such a 
basis in the general case. 

One of the main ideas of our approach is to construct a Riesz basis of finite- 
dimensional subspaces which are invariant for the generator of the semigroup 
(see (ED). 

In [3T| I32j we obtained the following general result. 

Theorem 3.2. There exists a sequence of invariant for A (see ^) finite- 
dimensional subspaces which constitute a Riesz basis in M-i- 

More precisely, these subspaces are {Vm \ \k\ > N,m — 1, ..,£} and a 2{N + 
l)n- dimensional subspace spanned by all eigen- and rootvectors, corresponding 
to all eigenvalues of A, which are outside of all circles Lm\ \k\ > N,m — 1, ..,£. 

Here vl''^ = P^^M2, where 

^-'^^^=2li4,^(-^'')'' 
are spectral projectors; Lm^ are circles defined before. 
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We emphasize that the operator A may not possess in a Riesz basis of 
generalized eigenspaces. We illustrate this on the following 

Example 3.3. Consider the particular case of the system |(T|): 



One can check that the characteristic equation is det A(A) — (a — A + Ae^'^)(/3 — 
A + Ae"'*') = and for a ^ (3 there are two sequences of eigenvectors, such that 
ll'^^n ~ ''^nll ~^ Oi as n ^ oo. It is clear that such family vectors do not form a 
Riesz basis. 



By stability we mean here asymptotic stability. For our neutral type system, 
as for several infinite dimensional systems, we have essentially two notions of 
asymptotic stability : exponential (or uniform) stability and strong stability. 

Definition 4.1. A linear system in a Banach space X is exponentially stable if 
the e-^* semigroup verifies: > 1, 3w > 0, Va;, ||e-^*a;|| < M;je~'^'[|a;||. 

The system is strongly stable if^x, Ije-^'a;!! ^0, as t ^ cg. 

The problem of exponential stability was widely described in several classical 
works. An sufficiently exhaustive analysis may be found in [39j (see also the 
references therein and the bibliographic notes). In our case the exponential 
stability is completely determinated by the spectrum of the operator A. It is a 
well known result for some linear neutral type systems: the spectrum has to be 
bounded away from the imaginary axis (cf. [TH Theorem 6.1]). 

Theorem 4.2. The system (0^ is exponentially stable if and only if a {A) C 
{A : ReA < -a < 0}. 

We can partially reformulate in terms of the matrix A-i the condition on 
the spectrum (t{A). 

Theorem 4.3. System (0) is exponentially stable if and only if the following 
conditions are verified 

i) a{A) C {A : RcA < 0} 



It can be interesting how the condition ii) of Theorem [331 may be formulated 
for the case of a general linear operator K in the system ^ . 

We would like to study more deeply the problem of asymptotic nonexponen- 
tial stability. To this end, we recall some important abstract result in this 
domain. We have the following 

Theorem 4.4. Let e^*, t>0 be aCo -semigroup in the Banach space X and A 
be the infinitesimal generator of the semigroup. Assume that {cr{A) O (iM.)) is at 
most countable and the operator A* has no pure imaginary eigenvalues. Then 
e-^* is strongly asymptotically stable (i.e. e-^*x — > 0, i — > +oo as x E X) if and 
only if one of the following conditions is valid: 





4 Stability 



li) a(A_i) C {A : |A| < 1}. 
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i) There exists a norm \\ ■ equivalent to the initial one \\ ■ \\, such that 
the semigroup e^* is contractive according to this norm: \\e"^'^x\\i < \\x\\i, 
yxG X,t> 0; 

a) The semigroup C^* is uniformly bounded: 3C > such that ||e'^*|j < C, 
t > 0. 

The Theorem 14.41 was obtained initially in [38] for a bounded operator A. 
The main idea were later used in [20] for the case of unbounded operator A, see 
also [2J for another approach. The proof in [20] follow the scheme of the first 
result [38]. The development of this theory concerns a large class of differential 
equations in Banach space (see [39] and references therein). A more genral 
result on the asymtotic behavior of the semigroup with respect to an arbitrary 
asymptotic was recently obtained in [35] . 

Our main result on the asymptotic stability of the neutral type system |(T])- 
([2]) is the following one. 

Theorem 4.5. Assume a{A) C {A : ReA < 0} and cti = a{A^i) n {A : |A| = 
1} 7^ Then the following three mutually exclusive possibilities exist: 

i) the part of the spectrum ai consists of simple eigenvalues only, i.e. to 
each eigenvalue corresponds a one- dimensional eigenspace and there are 
no rootvectors. In this case system is asymptotically stable. 

a) the matrix has a Jordan block, corresponding to ^ £ ai. In this case 
system ^ is unstable. 

Hi) there are no Jordan blocks, corresponding to eigenvalues in ai, but there 
exists /i £ (Ti whose eigenspace is at least two-dimensional. In this case 
system can be stable as well as unstable. Moreover, there exist two 
systems with the same spectrum, such that one of them is stable while the 
other one is unstable. 

The last case may be illustrated by a non trivial example (see also [32] for 
an example given partially in the M2-space framework). 

Example 4.6. (Rabah-Sklyar-Barkhaev 29 ) Consider the system 

z{t) - A^iz{t - 1) = Aoz{t) 

with 




= or 1. 



We have: <j{A) = {A : Ae'*' + A + e^ = 0} in C~, this can be proved by 
Pontriaguin Theorem [27]. The multipHcity of eigenvalues is clearly 2, and they 
do not depend of 7. The system is stable for 7 = and unstable for 7 7^ 0. 

5 Stabilizability 

We say that the system (2) is stabilizable if there exists a linear feedback control 
u{t) = F{zt{-)) — F{z{t + •)) such that the system (2) becomes asymptotically 
stable. 
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It is obvious that for linear systems in finite dimensional spaces the linearity 
of the feedback implies that the control is bounded in every neighbourhood of the 
origin. For infinite dimensional spaces the situation is much more complicated. 
The boundedness of the feedback law u — F{zt{-)) depends on the topology of 
the state space. 

When the asymptotic stabilizability is achieved by a feedback law which 
does not change the state space and is bounded with respect to the topology 
of the state space, then we call it regular asymptotic stabilizability. Under 
our assumption on the state space, namely 1, 0], C"), the natural Hnear 

feedback is 



Several authors (see for example [49^ \2E\ HU [13 and references therein) use 
feedback laws which for our system may take the form 



This feedback law is not bounded in 1, 0], C") and then stabiHzability 

is not regular. As a counterpart, they obtain exponentially stable closed loop 
systems. If the original system is not formally stable (see i-^- the pure 

neutral part (when A2 = A3 = 0) is not stable, the non regular feedback ([6]) 
is necessary to stabilize. From the operator point of view, the regular feedback 
law ([6]) means a perturbation of the infinitesimal generator A by the operator 
ej" which is relatively Abounded (cf. [TT]) and verifies V{A) = V{A + BT). 
Such a perturbation does not mean, in general, that A + BT is the infinitesimal 
generator of a Co-semigroup. However, in our case, this fact is verified directly 
[32^ [46j since after the feedback we get also a neutral type system like llj with 
V{A) = V{A + BT) (see below for more details). 

Prom a physical point of view, y^-boundedness of the stabiHzing feedback 
means that the energy added by the feedback remains uniformly bounded in 
every neighbourhood of (see also another point of view in [48]). Hence the 
problem of regular asymptotic stabilizability for the systems is to find a 

linear relatively yl-bounded feedback u = Tx such that the operator A + BT 
generates a Co-semigroup e(-^+^'^)* with V{A -f BT) = VIA) and for which 
||e('^+'^-^)*x|| ^ 0, as t ^ 00 for all x e 'D{A). The main contribution of this 
paper is that under some controllability conditions on the unstable poles of the 
system, we can assign arbitrarily the eigenvalues of the closed loop system into 
circles centered at the unstable eigenvalues of the operator A with radii such 
that X^^fe < This is, in some sense, a generalization of the classical pole 
assignment problem in finite dimensional space. Precisely we have the following 

Theorem 5.1. Consider the system (QJ) under the following assumptions: 

1) All the eigenvalues of the matrix A^i satisfy < 1. 

2) All the eigenvalues fij G cti (t(A_i) n {z : \z\ = 1} are simple (we denote 

their index j £ I). 




(6) 



where i^2(-),J^3(-) ei2(-l,0;C"). 



Y,F^z{t~h,)+ / F2{9)z{t + 9)dt+ / F3{9)z{t + 9)dt. 
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Then the system ^j) is regularly asymptotic stabilizable if 
3) rank (A^(A) = n for all RcA > 0, where 

A^(A) = -A/ + Ae-^A_i + A / e^M2(s)ds + / 



■0 



e^M3(s)ds, 



-1 



1 



4) rank - A_i B) = n for all \^\ = 1. 

6 Exact Controllability 

The problem of controllability for delay systems was considered by several au- 
thors in different framework. One approach is based on the analysis of time 
delay system in a module framework (space over ring, see [22]). In this case 
the controllability problem is considered in a formal way using different inter- 
pretations of the Kalman rank condition. Another approach is based on the 
analysis of time delay systems in vector spaces with finite or infinite dimension. 
A powerful tool is to consider a delay system as a system in a Banach functional 
space, this approach was developed widely in [13]. Because the state space for 
delay systems is a functional space, the most important notion is the function 
space controllability. A first important contribution in the characterization of 
null functional controllability was given by Olbrot [23] by using some finite 
dimensional tools as (A, i3)-invariant subspaces for an extended system. For 
retarded systems one can refer to [21] (and references therein) for the analysis 
of function space controllability in abstract Banach spaces. The case of neu- 
tral type systems with discrete delay was also considered in such a framework 
(see O'Connor and Tarn [25] and references therein). A general analysis of the 
time delay systems in infinite dimensional spaces is given in the book [6] where 
several methods and references are given. 

The problem considered in this paper is close to that studied in [25]. In this 
work the exact controllability problem was considered for neutral type systems 
with discrete delay using a semigroup approach in Sobolev spaces W^^'^ and a 
boundary control problem. 

We consider the problem of controllability for distributed delay system of 
neutral type in the space M2{—h, 0; C") = C" x L2{—h, 0; C") which is natural 
for control problems. 

The semigroup theory developed here is based on the Hilbert space model 
introduced in [8]. One of our result is a generalization of the result in |25j . 
The main non trivial precision is the time of controllability. We generalize the 
results given [16] for the case of a single input and one locaHzed delay (see also 
[HISl]). The approach developed here is different from that of [25j. Our main 
results are based on the characterization of controllability as a moment problem 
and using some recent results on the solvability of this problem (see [3] for the 
main tools used here). Using a precise Riesz basis in the space M2(— /i, 0;C") 
we can give a characterization of null-controllability and of the minimal time of 
controllability. 

The reachability set at time T is defined by 
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It is easy to show that TZti C TZt2 as Ti < T2. An important result is 
that TZt C 'D^A) C M2. This non-trivial fact permits to formulate the null- 
controllability problem in the following setting: 
i) To find maximal possible set TZt (depending on T)\ 

a) To find minimal T for which the set TZt becomes maximal possible, i.e. 
TZt=V{A). 

Definition 6.1. The system is said null- controllable at the time T if TZt = 

The main tool is to consider the null-controllability problem as a problem of 
moments. 



6.1 The moment problem 

In order to formulate the moment problem we need a Riesz basis in the Hilbert 
space M2. We recall that a Riesz basis is a basis which may be transformed 
to an orthogonal basis with respect to another equivalent scalar product. Each 
Riesz basis possesses a biorthogonal basis. Let {ip} be a Riesz basis in M2 
and {'ip} the corresponding biorthogonal basis. Then for each x G M2 we have 
X = YlipeM^^^^)^- ^ separable Hilbert space there always exists a Riesz 
basis. 

A state X = ^^q) ^ reachable at time T by a control u(-) G 

L2(0, T; C) iff the steering condition 



[ e^'Bu{t)dt. 
Jo 



(8) 



holds. This steering condition may be expanded using the basis {1^9}. A state x 
is reachable iff 

^ r {e^'Bu{t),^)dt^, 

for some u{-) G L2(— ft., 0; M''). Then the steering condition ^ can be substi- 
tuted by the following system of equalities 

T 

At 



{x,tP)= / {e-^'Bu{t),tP)dt, W- (9) 

"'0 

Let {61, . . . , br} be an arbitrary basis in luiB, the image of the matrix B and 
e M2, i = 1, . . . , r. Then the right hand side of ([9]) takes the form 

/ {e^'Bu{t),^P)dt = y2 [ (e-^*b„V}M,(i)di. (10) 
Jo ,,^1 Jo 

Effectiveness of the proposed approach becomes obvious if we assume that the 
operator A possess a Riesz basis of eigenvector. This situation is characteristic, 
for example, for control systems of hyperbolic type when A is skew-adjoint 
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{A* = —A) and has a compact resolvent (see, for example, [1], [16], [17]). Let 
in this case {tpk}, k £ N, he a orthonormal eigenbasis with Aipk — i^kVk, 
Xk € M. Assuming for simplicity r = 1, bi = b = ^f-O-k^k, oik 7^ 0, we have 
from (4), (5) 



Oik 



e 



u{t)dt, fc G N, (11) 



where x ~ J2k^kfk- Equalities (6) are a non- Fourier trigonometric moment 
problem whose solvability is closely connected with the property for the family 
of exponentials e~''^'=*, fc G N, to form a Riesz basis on the interval [0,T] ([1]). 
In particular, if e"''^''* forms a Riesz basis of L2[0, To] then one has 



TZt = {x :} [ — ] <oo} for all T > Tq. (12) 




Obviously formula (fT2l) gives the complete answer to the both items of the 
controllability problem. Returning now to neutral type systems we observe that 
the operator A given in is not skew-adjoint and, moreover, does not possess a 
basis even of generaHzed eigenvectors. So the choice of a proper Riesz basis to 
transform the steering condition in a moment problem is an essentially more 
complicated problem. 



6.2 The choice of basis 

In order to design the needed basis for our case we use spectral the properties of 
the operator A obtained in [32]. Let /ii, . . . , ^ he eigenvalues of A-i 

and let the integers Pm he defined as : dim (A_i — ^ml)" = Pm, m = 1, . . . 
Denote by 

A(^) \ (ln|^„| +i(arg^„ + 27rA:)) , m = 1,. fc G Z, 

n 

and let Lm"* be the circles of the fixed radius r < tq = | min | Am"* — x!f ^ \ centered 
at Am ^ . 

Let {Vm^} k£j^__ be a family of yl-invariant subspaces given by 

The following theorem plays an essential role in our approach 

Theorem 6.2. J31l\32^ There exists Nq large enough such that for any N > Nq 

i) dimVl''' = p„„ k>N, 

ii) the family \Vm ^} \k\>N U Vn forms a Riesz basis (of subspaces) in M2, 

m=i e 

where Vn is a finite- dimensional subspace (dim V^r — 2{N +l)n) spanned by all 
generalized eigenvectors corresponding to all eigenvalues of A located outside of 
all circles Lm\ \k\ > N, m — !,...,£. 

Using this theorem we construct a Riesz basis {ip} of the form 

{ift^^, |fc| > TV; m = 1, . . . , Z; J = 1, . . . ,pm} U {ipfj = 1,...,2{N+ l)n} 
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where for any m — and k : \k\ > N the collection {v?^ j}j=i,...,Pm in 

a special way chosen basis of Vm'^ and {^f}j=i,...,2{N+i)n is a basis of Vjv- In 
this basis equalities (4) with regard to (5) turns into a moment problem with 
respect to a special collection of quasipolynomials. Analyzing the mentioned 
moment problem by means of the methods given in [1] we obtain our main 
results concerning the null-controllability problem. 



7 The main results 

The characterization of the null-controllability is given by the following Theo- 
rem. 

Theorem 7.1. The system 0) is null-controllable by controls from L2(0,T) for 
some T > iff the following two conditions hold: 

i) rank[A^(A) B] — n, VA 6 C; where 

/Q pO 
-h J-h 

ii) rank [B A^iB ■ ■ ■ A'^Y^B] = n. 

The main results on the time of controllability are as follows. 

Theorem 7.2. Let the conditions i) and ii) of Theorem \7.1\ hold. Then 

i) The system is null- controllable at the time T as T > nh; 

ii) If the system 0) is of single control (r — 1), then the estimation of the 
time of controllability in i) is exact, i.e. the system is not controllable at 
time T = nh. 

For the multivariable case, the time depends on some controllability indices, 
suppose that dim_B — r. Let {&i, . . . , br} be an arbitrary basis noted f3. Let us 
introduce a set integers. We denote by Bi ~ . . . , 6^) , i = 0, 1, . . . , r — 1, 

which gives in particular Bq = B and i?r-i = br and we put formally Br = 0. 
Let us consider the integers 

nf = rank ■•■ A'Ti^B^.i], i = l,...,r, 

corresponding to the basis /3. We need in fact the integers 



Let us denote by 

B ■ B 

m-min = niax mV Wmax — nim max , 

P 13 i ' 

for all possible choice of a basis (3. 

The main result for the multivariable case is the following Theorem. 

Theorem 7.3. Let the conditions i) and ii) of the Theorem 1 7. 1\ hold, then 

i) The system ^ is null- controllable at the time T > m,„ax^; 
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ii) The system ^ is not null- controllable at the time T < nitninh. 

The proofs are based on the construction of a special Riesz basis of A- 
invariant subspaces in the space M2 according to [3T| [32] and on the analysis of 
the properties of some quasi-exponential functions to be a Riesz basis in ^2(0, T) 
depending of the time T [3]- 
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